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ABSTRACT 


The  achievable  gain  in  the  radiation  resistance  and  directivity  of  a 
low  frequency  underwater  transducer  due  to  the  presence  of  an  array 
of  sympathetic  resonators  was  analyzed.  The  resonators  were  all  taken 
to  be  air  bubbles,  and  both  the  resonators  and  transducer  were  taken 
to  be  compact  (ka«l).  The  resonators  were  taken  to  be  equally 
spaced  around  a  circle  of  radius  R,  with  the  transducer  located  on  the 
axis.  The  gain  was  calculated  for  various  numbers  of  resonators  as  a 
function  of  kare9onator.  katretnsducer  and  kR,  for  the  transducer  in  the 
plane  of  the  resonators  and  out  of  the  plane  a  distance  of  one-quarter 
wavelength.  For  the  transducer  in  the  plane,  a  gain  in  radiation 
resistance  of  approximately  two  is  possible  with  six  or  more 
resonators.  For  the  transducer  out  of  the  plane,  it  is  shown  that  a 
significant  gain  in  directivity  can  be  achieved  at  the  expense  of  a  slight 
decrease  in  the  gain  in  the  radiation  resistance. 
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I.  INTRODUCTION 


A.  BACKGROUND  AND  HISTORY 

I.  Tolstoy  [Ref.  1], [Ref  2]  and  A.  Tolstoy  and  I.  Tolstoy  [Ref.  3]  nave 
recently  analyzed  the  multiple  scattering  of  sound  from  a  collection  of 
compact  (ka«l)  monopole  scatterers  (e.g.  gas  filled  bubbles,  gas  filled 
balloons  and  thin  shells  in  water)  insonified  by  a  simple  harmonic 
plane  wave  at  frequencies  close  to  the  resonance  frequency  of  an 
individual  scatterer.  The  response  of  these  systems  was  quantified  by 
computing  the  ratio  of  the  scattered  pressure  on  the  surface  of  a 
scatterer  in  the  presence  of  the  other  elements  to  that  of  the  incident 
plane  wave.  This  ratio  was  defined  as  the  amplification  factor,  |i.  They 
investigated  the  value  of  this  amplification  factor  for  various  simple 
geometrical  arrangements  of  monopole  resonant  scatterers,  and 
showed  that  amplification  factors  of  the  order  of  102  were  to  be 
expected  for  air  bubbles  in  water  at  one  atmosphere  ambient  pressure. 
(The  amplification  factor  of  a  single  bubble  is  about  70.)  These 
systems  were  said  to  be  "quasiresonant". 

I.  Tolstoy  [Ref.  4]  has  also  recently  analyzed  the  strength  and 
directionality  of  the  resonant  multiple  scattering  of  sound  from  a 
linear  array  of  compact  monopole  resonant  scatterers  irradiated  by  a 
plane  wave.  He  reports  that  "quasiresonance"  effects  exhibit  strong 
directionality  in  this  system,  with  amplification  factors  of  up  to  700. 
The  potential  application  of  resonant  multiple  scattering,  or 
"quasiresonance",  in  Tolstoy's  words,  to  improve  the  performance  of  a 
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low  frequency  active  sonar  transducer  is  the  subject  of  the  research 
described  in  this  thesis. 

B.  OBJECTIVES 

The  objectives  of  this  research  are  multiple.  The  first  is  to 
reproduce  I.  Tolstoy’s  results  using  network  analysis,  a  tool  familiar  to 
the  transducer  and  array  designer.  It  will  then  be  shown  that  the  two 
methods  of  analysis,  multiple  scattering  and  network  analysis,  give 
identical  results.  The  existence  of  quasiresonance  will  then  applied  to 
a  circular  array  of  scatterers,  where  the  incident  plane  wave  is 
replaced  by  an  active  transducer  located  on  the  axis  of  the  circle  and 
radiating  spherical  waves.  The  acoustic  advantage  gained  in  this 
system  will  be  quantified  by  the  gain  in  the  radiation  resistance  seen 
by  the  transducer  in  the  presence  of  these  scatterers.  An  increase  in 
radiation  resistance  enables  the  transfer  of  more  power  to  the  acoustic 
field  for  a  given  surface  velocity  of  the  transducer.  (High  power  low 
frequency  active  sonar  transducers  tend  to  be  resonant  devices  that 
are  displacement  limited.)  The  cases  of  the  transducer  in  the  plane  of 
the  circle  and  displaced  by  one-quarter  wavelength  from  the  plane  will 
be  analyzed  for  both  their  increase  in  radiation  resistance  and  gain  in 
directivity  over  that  of  a  single  transducer. 
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II.  THEORY 


A.  QUASIRESONANCE 

For  simple  configurations  of  compact  resonant  scatterers 
insonified  by  a  simple  harmonic  plane  wave,  Tolstoy  predicted  the 
existence  of  so-called  "quasiresonance"  using  multiple  scattering 
theory.  (Quasiresonance  is  Tolstoy's  term  for  the  enhancement  of  the 
scattered  pressure  from  a  resonant  scatterer  as  a  result  of  multiple 
scattering.)  Tolstoy's  analysis  will  be  reviewed  and  applied  to  the 
simplest  case  in  section  Al.  In  section  A2  the  same  problem  is  treated 
using  network  analysis.  The  existence  of  quasiresonance  is  again 
predicted  using  a  coupled  network  model.  It  will  be  shown  that  this 
second  method  produces  results  identical  to  Tolstoy's. 

1.  TOLSTOYS  ANALYSIS  OF  MULTIPLE  SCATTER 

Tolstoy's  analysis  of  the  multiple  scattering  of  a  plane  wave 
incident  upon  a  collection  of  compact  resonant  scatterers  proceeds 
using  linear  theory  of  acoustics.  The  following  development  parallels 
that  of  Reference  [1]  and  Reference  [3],  with  the  difference  being  that 
here  all  time-harmonic  functions  are  taken  to  be  of  the  form  elwt.  This 
will  result  in  minor  differences  in  the  appearance  of  otherwise 
equivalent  equations. 

Consider  first  the  case  of  a  single  resonant  scatterer  insonified 
by  an  incident  plane  wave.  From  Lighthill  [Ref.  5),  the  sound  field  of  a 
spherical  source  is  given  by  Equation  2-1. 
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p(r,t)=  (47tr }" 1  q(t-r/c) 


(2-1) 


where: 

r  =  Distance  from  scatterer  center  to  the  field  point, 
q  =  Mass  injection  rate  of  a  spherical  sound  source, 
c  =  Speed  of  sound  in  water. 


The  dot  over  q  indicates  the  derivative  with  respect  to  time.  Treating 
the  single  scatterer  as  an  air  fillec  Dubble,  for  an  incident  plane  wave 
of  free  field  pressure  pf,  q  takes  the  form  of  Equation  2-2. 


q  =pfk:V 


ypp  { | 

|l  -U-jka)e4| 

pc2  ! 

1  «i' 

(2-2) 


where: 

pf  =  Incident  free  field  acoustic  pressure, 
k  =  Wave  number,  oVc. 

V  =  Volume  of  scatterer. 
a  =  Scatterer  radius, 
y  =  Ratio  of  specific  heats  (air), 
c  =  Speed  of  sound  in  water, 
p  =  Density  of  water. 

P0  =  Ambient  pressure  of  environment. 

(o  =  Angular  frequency  of  plane  wave. 
coq  =  Resonance  frequency  of  scatterer. 


Inserting  Equation  2-2  into  Equation  2-1,  dropping  the  time 
dependence,  and  solving  for  the  ratio  of  the  scattered  acoustic 
pressure  at  some  radius  r  from  the  scatterer,  pfingie(r),  to  that  of  the 
incident  pressure,  pf  ,  one  obtains  Equation  2-3. 


Psinglc(r)  _  J  |J-.j  e-jkr 

pf  Mi-jd 
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(2-3) 


where: 


(to^/to2  -l) 
ka 


For  situations  where  the  surface  tension  and  corrections  due  to  energy 
losses  in  the  interior  of  the  scatterer  (e.g.  bubbles,  balloons)  are 
negligible,  the  ratio  of  the  resonance  angular  frequency  of  the 
scatterer  to  the  angular  frequency  of  the  incident  plane  wave  is  given 
by  Equation  2-4. 


ojo  _  /  3yPo 

03  A/  pc2(kap 


(2-4) 


Letting  r  equal  a,  the  scatterer  radius  in  Equation  2-3,  results  in 
Equation  2-5,  the  ratio  of  the  scattered  acoustic  pressure  at  the  single 
scatterer's  surface  to  that  of  the  incident  free  field  pressure. 


nsc 

Fsingle  -  _L  l_L_\  e-.jka 
pf 


(2-5) 


Tolstoy  introduces  the  quantities  Fa  and  A  by  rewriting  Equation 
2-3  as  Equation  2-6. 


where: 


A  =  -i 


(2-6) 


5 


Fa  is  the  normalized  scattered  acoustic  pressure  for  a  single  scatterer 
insonified  by  an  incident  plane  wave  and  A  is  the  single-scatterer 
amplitude  coefficient. 

Now  consider  a  system  of  N  identical  interacting  scatterers 
insonified  by  a  plane  wave.  Let  p^c(r)  represent  the  scattered  acoustic 
pressure.  An  equation  for  p£c(r)/pn,  equivalent  to  that  of  Equation  2-6 
for  a  single  scatterer,  can  be  written  as  Equation  2-7  for  the  multiple 
scatterer  situation. 


rb  =e0iUb„^ 


Pn 


jkr 


(2-7) 


rbn  is  the  normalized  scattered  acoustic  pressure  from  scatterer  n, 

insonified  by  both  the  incident  plane  wave  and  the  fields  of  the  other 
scatterers,  Pn  is  the  free  field  pressure  of  the  incident  plane  wave  at 
scatterer  n,  and  Bn  is  the  multiple  scatterer  amplitude  coefficient  for 
scatterer  n. 

For  compact  scatterers,  the  overall  response  of  a  scatterer  to 
multiple  incident  pressure  fields  can  be  evaluated  by  summing  up  the 
response  to  each  individual  incident  pressure  field  (Ref.  3].  For  a 
system  of  N  scatterers,  this  results  in  Equation  2-8  for  the  scattered 
field  at  the  surface  of  the  n*h  resonator. 


B  _  a 

Dn  i  — 

jkan 


1  +  1  Br 


,-jkJr, 


jkl„ 


e'N'n. 


m*n 


g-jkan 

jkan 


(2-8) 


where:  is  the  phase  difference  in  the  incident  plane  wave  between 

resonator  m  and  n,  and  m=l...N.  The  subscript  notation  for  A  has 
been  added  to  allow  for  a  system  of  different  sized  scatterers. 


Tolsto>  defines  the  amplification  factor  p  (different  than  the 
scatter  amplitude  coefficient  mentioned  previously)  to  be  the 
magnitude  of  the  ratio  of  the  scattered  acoustic  pressure  on  the 
surface  of  a  scatterer,  PnC,  to  the  incident  free  field  acoustic  pressure 
of  the  incident  plane  wave,  p„.  With  this  definition,  Equation  2-9 
follows  from  Equation  2-7. 


Mn  = 


IPtl 

Pn 


|Bnj 

kan 


(2-9) 


Tolstoy  denotes  the  amplification  factor  for  a  single  ocatterer  as  g0. 
Equation  2-10,  then,  follows  from  Equation  2-6. 


Mo 


n 


Pn 


kan 


(2-10) 


From  Equations  2-9  and  2-10,  g/po  is  expressed  in  terms  of 
Bn/An  by  Equation  2-11. 


Mn 


Mo 


n 


=  |Bn/Anl 


(2-11) 


a.  APPLICATION  OF  TOLSTOY’S  METHOD 

As  an  example  application,  Tolstoy  considered  a  system  of  two 
identical  resonators  insonified  by  a  plane  wave  traveling  along  the  axis 
joining  their  centers,  as  shown  in  Figure  2-1. 
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Resonator  of 


Figure  2-1.  A  Two  Resonator  System  Insonified  by 
a  Plane  Wave  Along  the  Axis  of  The  Resonators. 


Application  of  Equation  2-8  to  the  situation  shown  in  Figure  2-1 
results  in  the  set  of  Equations  2-12. 


Bj  e'jkn  =A^-  +  A  B-,  ^  e-i« 
jka  jka  jka  *  jkl 


B^=A^  +  A^B1^eW 

jka  jka  jka  jkl  (2-12) 


For  a  system  of  identical  scatterers,  the  values  of  A!  and  A2  are  equal. 
Therefore,  the  subscript  on  A  in  Equation  2-12  has  been  omitted. 
Since  all  variables  in  Equation  2-12  are  known  except  Bj  and  B2,  one 
can  solve  for  Bi/A  in  terms  of  A  and  kl,  resulting  in  Equation  2-13. 


B,/A  = 


1  + 


e-jkJ 


(2-13) 


Substituting  the  expressions  for  £  and  to/ coo  for  an  air-filled  bubble  into 
the  expression  for  A,  B\/A  is  given  in  terms  of  the  variables  ka,  kl,  c, 
Po,  and  p  by  Equation  2-14. 


B,/A  = 


1  -  e  Jkl 

kaej(ka-ki)  \ 

kl 

oc  ((ka^+ilka 

| 

Pf  _ 

1 

P  pc2(ka) 

II 

nsc 

Fsinglc 

l  - 

1  kiLgj(ka-ld)  \" 

'  kl 

> 

ipc 

(ikap+jlka  -  ^  )| 

1  1  pc^ka)  II  1 

(2-14) 


For  bubbles  just  under  the  surface  of  the  ocean,  we  take: 

P0  =  10132  Pa  (1  atmosphere) 

p  =  998  kgnr3 
c  =  1498  ms1 

As  observed  by  Tolstoy  (Ref.  3J,  plots  of  Equation  2-14  display 
spurious  peaks.  Following  Tolstoy,  then,  the  graphs  of  IBi/AI  versus 
ka  and  kl,  plotted  in  Figures  2-2a,b  have  been  treated  to  limit  peaks  in 
bound  intervals  less  than  2ka.  (According  to  Tolstoy  peaks  narrower 
than  2ka,  where  2a  equals  the  resonator  diameter,  are  not  meaningful, 
since  the  pressure  over  a  compact  resonator  is  approximately 
constant.)  The  truncation  method  used  involved  comparing  three 
consecutive  values  of  IBi/AI  in  the  kl  direction  for  each  value  of  ka. 
(Three  consecutive  points  in  the  kl-direction  span  approximately 
2ka.)  If  the  value  of  IBi/AI  for  the  middle  point  was  greater  than 
either  of  its  neighboring  points  (a  peak)  then  this  value  was  set  equal 
to  the  mean  of  the  two  neighboring  points.  This  procedure  was 
repeated  several  times  to  get  the  desired  smoothing. 
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Figures  2-2a,b  show  a  maximum  response  of  IBj/AI  equal  to  5.5 
at  approximately  ka=0. 01393  and  kl=0.60.  As  a  result  of  multiple 
scattering,  then,  the  scattered  pressure  amplitude  of  the  first  bubble 
is  5.5  times  greater  in  the  presence  of  the  second  bubble  than  in  its 
absence.  Additionally,  the  value  of  ka  at  the  maximum  response  is 
shifted  from  that  for  a  single  resonant  bubble  of  .01379. 


2.  NETWORK  ANALYSIS  OF  MULTIPLE  SCATTER 

Fig.  2-3  displays  a  network  model  of  N  interacting  scatterers 
insonified  by  an  incident  acoustic  wave.  (Note:  pressure  (p)  and 

volume  velocity  (Un)  are  chosen  as  the  mechanical  variables,  and  the 
direction  of  positive  velocity  is  clockwise.)  The  network  equations  for 
Fig.  2-3  are: 

N 

Pn  =  Pn  +  ^  Z‘nrmUm  P„  =  -Z'n  Un  (2-15) 

ni  =  I 


where  for  Fig  2-3: 

Pn=  Surface  pressure  at  element  n. 

p/,=  Free  field  pressure  at  element  n  due  to 
an  externally  incident  field  only. 

Un=  Volume  velocity  of  element  n. 

Z^=  Acoustic  impedance  of  element  n. 

Z^n=  Acoustic  radiation  impedance  seen  by 
element  n  with  all  other  elements 
blocked.  (Um*n=0  m=  1  ...N) 

7^n=  Acoustic  transfer  impedance  from 
element  m  to  element  n. 


Resonators 


Acoustic  Field 


Figure  2-3.  Network  Model  of  N  Interacting 
Resonators  Insonified  bv  a  Plane  wave. 


Neglecting  the  effects  of  surface  tension,  thermal  relaxation, 
etc.,  the  intrinsic  acoustic  impedance  of  a  small  (ka«l)  gas-filled 
bubble  can  be  approximated  by  one  of  reactance  only,  given  by  Kinsler, 
Frey,  Coppens  and  Sanders  IRef.  6]  as  Equation  2-16. 


where: 


Zan=j 


-3ypn  s 
dkaj 


- 1 
n 


an  =  Radius  of  element  n. 
c  =  Speed  of  sound  in  water, 
k  =  Wave  number,  oVc. 

P(>  =  Ambient  pressure  of  environment. 


(2-16) 
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Sn  =  Surface  area  of  element  n. 

7  =  Ratio  of  specific  heats  (air). 

The  self  acoustic  radiation  impedance  seen  by  element  n,  Z‘nn. 
can  be  approximated  by  Equation  2-17. 

Z|in  =  PC  (|k£ln)  Sn  (2-17) 

where: 

P  =  density'  of  water. 

Equation  2-17  is  the  acoustic  radiation  impedance  presented  to  a 
compact  source  with  no  other  active  or  scattering  elements  present. 
It  is  also  the  acoustic  radiation  impedance  seen  bv  a  compact  source  in 

the  presence  of  other  compact  sources  when  these  sources  are 
blocked  (Un1;eil  =  0  m=l...N). 

The  acoustic  transfer  impedance  from  element  m  to  element  n 
is  given  by  Equation  2-18. 


Zar  _  7. u 

11m  “  ^mm 


kiln. 


kl„ 


ej(ka,n-kl„.„) 


(2-18) 


In  matrix  form  the  network  equations  (Equation  2-15)  can  be 
written  as  Equation  2-19. 


!pni=  teJjUml+lpj,!  m.n  =  1...N  (2-19) 

where: 

| Xn  |  n=  1...N  A  lxN  column  matrix,  N  elements. 

'  X  n  ni  n.m  =  I...N  A  NxN  square  matrix.  N2  elements. 

Substituting  pn  =  -Z^U,,  and  rearranging  results  in  Equation  2-20. 
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where 


!-pU=  [Znm+Z^5nm]  { Um }  m,n  =  1...N 


(2-20) 


Snm  =  Kronecker  delta  function, 
|  1;  m=n  I 

5nm  = 

|  0;  m*n  J 


a.  APPLICATION  OF  NETWORK  ANALYSIS 

Modeling  the  system  of  two  scatterers  shown  in  Figure  2-1  as  a 
network.  Equation  2-20  becomes  Equation  2-21. 


-Pi 


-P2 


Zar  ,  ^7 a 

Zar  ^ar  ,7a 

2 1  ^22  '  ^2 


u, 

u2 


(2-21) 


Using  the  relation  p^c  =  ZnnUn  to  eliminate  Un,  the  above  equations 
can  be  written  in  terms  of  the  scattered  acoustic  pressure  at  each 
individual  resonator,  pnC,  and  the  incident  plane  wave  pressures  as 
given  bv  Equation  2-22. 

P*c  =  (pfi  +  C12  p|c  )  Dj1 

PS2C  =  (pf2  +  C21  pf  )  D2‘  (2-22) 


where 


Cnm  — 


yar 

^nm 

SSn, 


Dn  = 


Zanrn  +Z* 

Zrin 


For  a  system  of  identical  resonators,  Cnm  and  Dn  are  the  same  for  all  m 
and  n;  therefore  the  subscript  notation  is  not  required.  Solving 
Equation  2-22  for  the  scattered  acoustic  pressure  at  each  resonator  in 
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terms  of  the  incident  plane  wave  pressure  results  in  Equation  2-23 
for  the  scattered  pressure  at  resonator  one. 


,c  P^P^C 

1  D2  -  C2  (2-23) 

The  scattered  acoustic  pressure  at  the  first  resonator  in  the 
absence  of  the  second  resonator,  plingiep  is  obtained  using  Equation  2- 

23  by  letting  the  transfer  coefficient  C  equal  zero,  resulting  in 
Equation  2-24. 

nsc  =  R1 

Psmg'cj  D  (2-24) 


Dividing  equation  Equation  2-23  by  2-24  results  in  Equation  2- 
25  for  the  ratio  Bi/A  as  defined  by  Tolstoy. 


B  i/A  = 


(2-25) 


Inserting  Equations  2-16  thru  2-18  into  Equation  2-25,  one 
obtains  an  equation  for  Bi/A  in  terms  of  the  variables  ka,  kl.  c,  Po.  pfi  . 
P2  and  P  as  given  by  Equation  2-26. 
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i  -  (pS/pi) 


B!/A  =  -  Pj--  = 

nsc 

V  single] 


kSLgjlka-kl) 

M _ 


|  pc  |(kap  +j| 


ka  - 


3yPq 

pc2(ka) 


1  - 


kg-jlka-kl) 

kT 


(2-26) 


pc 


j(ka)2  +j| 


ka  - 


3yPo 


pc^ka)  II  I 


f  f 

The  incident  plane  wave  pressures  pi  and  P2  are  of  equal  magnitude 
and  differ  only  by  the  phase  factor  el^12,  resulting  in  Equation  2-27. 


P2  =  e-jkl 


Pi 


(2-27) 


Substituting  this  expression  into  Equation  2-2G,  one  obtains  Equation 
2-28. 


Eh/A  = 


1  -e-iw 


kjM(ka-kl) 

kl 


1  - 


ka^i(ka-kl) 
Jd _ 


|  pc  J(kap  +jj 


ka 


3yP0  1 

pc^ka)  if  I 


(2-28) 


Equation  2-28  is  identical  to  Equation  2-14,  demonstrating  the 
equivalence  of  Tolstoy's  method  and  the  network  analysis  method. 


3.  A  TRANSDUCER  IN  A  MULTI-RESONATOR  FIELD 

Network  analysis  will  now  be  applied  to  the  problem  of  a  compact 
transducer  in  the  presence  of  an  array  of  interacting,  compact 
resonant  scatterers.  The  procedure  is  the  same  here  as  in  section  A2, 
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except  that  the  incident  plane  wave  is  replaced  by  an  outward 
traveling  spherical  wave  whose  source  is  a  compact  transducer.  A 
network  model  of  the  transducer  and  N  resonator  system  is  shown  in 
Figure  2-4. 


Acoustic  Field 


Transducer 


Resonators 


Fig.  2-4.  Network  Model  of  a  Single  Transducer 
and  an  N  Resonator  Field. 


The  network  equations  for  Fig.  2-4  are; 

Pn  =  X  Z^niUm  pn  =  -Z„  Un  for  n  =  1  ...N  Po  =  TI  -  Zao  U()  (2-29) 

m 


where 

T  =  transduction  coefficient. 

Zao=  transducer  open  circuit  (1=0)  acoustic  impedance. 
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and  the  subscript  n=0  denotes  the  transducer,  n*0  a  resonator. 

The  transducer's  open  circuit  acoustic  impedance  (Zao)  and 
transduction  coefficient  (T)  are  shown  only  for  completeness  and  will 
not  be  considered  further,  as  the  quantity  of  interest  will  be  the 
radiation  impedance  presented  to  the  transducer. 

Using  the  network  Equations  2-29,  the  matrix  equation  for  the 
network  model  shown  in  Figure  2-4  is  given  by  Equation  2-30. 

Po 

<n  ' 
iPn  / 

where 

!xno)  n=l...N  denotes  a  lxN  column  matrix,  N  elements. 

[Xnm]  n,m=l...N  denotes  a  NxN  square  matrix,  N2  elements. 

and  the  superscript  T  denotes  transpose.  Replacing  the  1  x  N  surface 

pressure  matrix,  !pn),  with  l-ZnUn}  and  then  normalizing  by  U0,  the 

transducer's  volume  velocity,  Equation  2-30  becomes  Equation  2-31. 

yar 
^total 

in' 

»u/ 

where 

ztotai  =  Total  acoustic  radiation  impedance  seen  by  the 
transducer. 

Um1  =  Ratio  of  element  m  s  volume  velocity  to  that  of  the 
transducer  (Um1  =  Um/U«) . 

Equation  2-31  can  be  written  as  two  separate  matrix  equations, 
Equations  2-32a,b. 


yar 

MKI 

/yar  |T 
\AnO| 

1 

iy  ar  1 

l^nOf 

[Z%m  +ZS5nm] 

lUR1! 

;  n,  m  =  1...N 


(2-31) 


7ar  1 7 ar  |T 


-00 

/yar  | 
lAiOl 


\  m()  | 


[Z; 


ar  | 
ran: 


U0 

{ Um } 


n,  m  =  1...N 


(2-30) 
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Z“,»i  =  Zoo+'Zrf!T|USl]  (2-32a) 

{0|  =  |z"„]  +  [zss,,  +zj6„m](u;;'!  (2.32b) 

Solving  Equation  2-32b  for  fUm1}  by  matrix  inversion  and  substituting 
the  result  into  Equation  2-32a,  one  obtains  Equation  2-33,  the  total 
acoustic  radiation  impedance  seen  by  the  transducer. 

zf„„,  =  ZR|  ■  Ko!T  [zj.  +Z;8„J-‘  |Z"n'„!  (2-33) 

Equation  2-33  expresses  this  impedance  in  terms  cf  the  radiation  and 
transfer  impedances  of  the  individual  elements. 

The  tot'd  acoustic  power  radiated  by  the  transducer  is  given  by 
Equation  2-34. 


n  =  i|U,J2Re(Z"„l 


(2-34) 


The  gain  in  power  delivered  to  the  acoustic  field  due  to  the  addition  of 
the  sympathetic  resonators  equals  the  ratio  of  total  power  delivered 
with  the  resonators  present  to  the  power  delivered  when  they  are 
absent,  n  =  ~\  U<J"  Re(Z("u  ).  For  a  displacement-limited  source,  typical  of  a 
sonar  projector  operated  near  its  resonant  frequency,  this  ratio 
reduces  to  the  ratio  of  the  radiation  resistances  for  the  two  cases.  This 
gain  factor  is  expressed  as  Equation  2-35. 

Gr  = 

Re(Z?x> )  (2-35) 
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m  .  ANALYSIS 


A.  ANALYSIS  OF  RESONATOR  AND  TRANSDUCER  CONFIGURATIONS 

The  theory  developed  in  the  previous  chapter  will  now  be  applied 
to  two  specific  geometric  configurations  of  a  transducer  and  N 
identical  resonators.  The  resonators  will  be  treated  as  simple  gas  filled 
bubbles.  It  is  assumed  throughout  that  the  acoustic  wavelength  is 
much  larger  than  the  dimensions  of  the  transducer  and  the 
resonators.  For  initial  calculations,  ambient  pressure  is  taken  to  be  one 
atmosphere,  and  the  radius  of  the  transducer  is  taken  to  be  equal  to 
that  of  the  resonators  (kao  =  kan)- 

1.  PLANAR  CONFIGURATION 

The  geometry  of  the  first  arrangement  to  be  analyzed  is  shown 
in  Figure  3-1.  It  consists  of  N  identical  resonators  equally  spaced 


I 

I*'  1 


Figure  3-1.  Resonators  Equally  Spaced  around  Transducer 
in  a  Planar  Configuration  of  Radius  lno.  N=4. 
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around  a  circle  of  radius  ln0.  This  geometric  arrangement  will  be 
referred  to  as  the  planar  configuration. 

For  the  configuration  shown  in  Figure  3-1,  Equation  2-35 
becomes  Equation  3-1. 

c  -  Re(  -  iZSolT  [ZS.  -|-Z5Snm1"'  |zy„j)  nm=)  4 

Re(  ZS,I  '  (3-1) 

Using  Equations  2-16,  17,  and  18,  and  canceling  the  common  surface 
areas  Sn.  Equation  3-1  can  be  written  as  Equation  3-2  in  terms  of 
specific  acoustic  impedances: 


Gr  = 


Re(  ?! 


HK) 


-  !z 


ar  lT 
mo  ( 


lznni  +zn5nm ' 


/_ar  \ 
<zn()> 


Re(  Z?)( 


)()’ 


n,m=i  ...4 


13-2) 


where: 


Z(H)  =  pc  ({ka|  --rjk.a! 

znm  +z?,Snn1  =  pc  ((ka)2+jka)  -  j  when  n=m 

c(ka) 

znm  +znSnm  =  znm  =  pc  ((ka)2+jka)  rj^-~  ej  (ka-klpj^)  when  n*m 

klrun 

z-;r„  =  pc  ((ka)2+jka!  ej  (ka-klno) 

kln0 

and  lnm  is  the  distance  between  element  n  and  m,  the  subscript  0 
denoting  the  transducer.  Using  simple  trigonometry,  klnm  for  n.m  *  0 
can  be  written  in  terms  of  m,  n,  and  klno  as  Equation  3-3  for  the 
planar  case. 


klnm  =  kln02  sin(^-|n-mj| 


where  n.m  =  1  ...N 


(3-3) 
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Using  Equations  3-2  and  3-3,  the  gain  in  radiation  resistance  seen  by 
the  transducer,  Gr,  can  be  expressed  in  terms  of  ka  and  kln0. 

The  gain  in  radiation  resistance  for  the  N=4,  planar 
configuration  can  be  plotted  in  the  form  of  surface  and  contour  plots 
using  Equation  3-1,  as  shown  in  Figures  3-2a  and  2b.  It  should  be 
pointed  out  that  here  and  in  all  subsequent  surface  plots  of  Gr,  no 
spurious  peaks  were  observed,  and  no  smoothing  has  been  applied. 
Inspection  of  these  graphs  reveals  two  local  maxima,  one  of  magnitude 
4.34  at  ka=0.975  kar,  klno=0.5  and  another  of  magnitude  1.79  at 
ka=  1.005  kar,  kln0=2.7,  where  kar  equals  0.01379,  the  resonance 
value  of  ka  for  a  single  air  filled  bubble  at  one  atmosphere  ambient 
pressure.  The  local  maximum  at  ka=0.975  kar,  kln0=0.5  is  an  artifact 
of  the  range  of  kln0  chosen.  In  fact,  the  gain  in  radiation 
resistance  diverges  as  (kln0)  1  when  kIno  approaches  ^ero.  This  case 
is  uninteresting  and  will  not  be  considered  further.  The  interesting 
local  maximum  at  ka=  1.005  kar,  klno=2.7  represents  a  gain  in  the 
radiation  resistance  seen  by  the  transducer  by  a  factor  of  almost  two. 
Note  that  kl  =  7t  corresponds  to  a  separation  distance  of  one-half 
wavelength.  A  maximum  gain  in  radiation  resistance  is  therefore 
obtained  when  the  resonators  are  located  a  distance  just  under  one- 
half  wavelength  from  the  transducer. 

The  gain  in  radiation  resistance  for  the  planar  configuration 
where  N  is  set  equal  to  eight  and  fifteen  is  shown  in  Figures  3-2c  thru 
3-2f.  These  plots  display  similar  characteristics  to  the  N  =  4  case.  All 
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show  a  local  maximum  of  magnitude  approximately  two  at  ka  =  kar  = 
0.01379  and  klno  =  x. 

The  local  maximum  values  for  Gr  at  kln0  =  k,  ka  =  kar  = 
0.01379  for  various  values  of  N  ranging  from  2  to  25  and  their 
associated  parameters  are  tabulated  in  Table  3-1.  Two  major  points 
can  be  observed  from  Table  3-1.  First,  the  resonant  bubbles’  volume 
velocity  leads  the  transducer's  by  a  phase  shift  of  approximately  90° 
(equivalent  to  -270°).  This  phase  shift  is  created  in  two  parts.  Since 
the  separation  between  transducer  and  resonator  is  approximately  one 
half  a  wavelength,  this  introduces  a  shift  of  -180°  in  phase  due  to 
simple  wave  propagation.  The  remaining  -90°  is  introduced  at  the 
surface  of  a  compact  resonating  bubble  between  the  incident  and 
scattered  wave  as  predicted  by  Equation  2-24.  Second,  note  that  little 
is  gained  as  one  increases  the  number  of  resonators  beyond  about  four. 


N 

Gr 

k^max 

klmax 

mag(U{f‘) 

argdJJf1) 

2 

Bia 

MM 

2.5 

.358 

83° 

4 

pH 

Pi  111 

2.7 

.546 

77° 

6 

11 

2.8 

.472 

92° 

8 

19 

imwm 

2.9 

.352 

88° 

10 

Pll 

.01386 

2.9 

.284 

90° 

15 

■  I 

.01386 

2.8 

.191 

94° 

20 

mm 

.01379 

2.9 

.144 

93° 

25 

m 

.01379 

2.8 

.113 

91° 

Table  3-1.  Values  of  Gr  for  Values  of  N  Ranging  from  N=2 
to  25  and  the  Parameters  Associated  with  the  Local 
Maximum  at  kln0  =  n,  ka  =  kar=0. 01379  for  the  Planar 
Configuration. 
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a.  PLANAR  CONFIGURATION  BEAM  PATTERN 

Polar  plots  of  the  beam  patterns  can  be  developed  for  the  planar 
configuration  using  the  axis  system  defined  in  Figure  3-3,  The 
origin  is  located  at  the  transducer,  with  the  resonators  located  in  the 
x,v  plane  and  0  being  the  angle  measured  from  the  positive  y  axis  in 
the  x,v  plane. 

Figure  3-4  clearly  demonstrates  the  acoustic  advantage  of 
the  planar  configuration  over  that  of  a  lone  transducer.  The  zero  dB 
reference  is  taken  as  the  field  of  an  identical  transducer  with  equal 
volume  velocity  in  absence  of  the  resonators. 


Figure  3-3.  Axis  System  for  Beam  Patterns  Plots,  Planar 
Configuration.  Note  the  definitions  of  v,z  and  z,  0=45°  Plane 
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z  axis 


Figure  3-4.  Beam  Pattern  for  the  Planar 
Configuration  in  the  y.z  Plane  and  the  z,0=45° 

Plane.  P0=  latm.  Scale  in  Decibels  (dB). 

2.  CONICAL  CONFIGURATION 

The  second  geometric  arrangement  analyzed  was  chosen  to 
take  advantage  of  the  90°  phase  shift  between  the  volume  velocity  of 
the  resonators  and  the  transducer  for  maximum  radiated  power.  If  the 
transducer  is  offset  from  the  plane  of  the  resonators  by  one  quarter 
wavelength  along  the  positive  z  axis,  the  radiation  from  the  resonators 
will  constructively  interfere  with  that  of  the  transducer  along  the 
positive  z  axis.  This  configuration  will  increase  the  directivity  of  the 
system.  This  arrangement  will  be  referred  to  as  the  conical 
configuration  and  is  shown  in  Fig.  3-5. 
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Figure  3-5.  Resonators  Equally  Spaced  Around  Transducer 
in  a  Conical  Configuration  at  Radius  R  and  N  =  4.  The 
Transducer  is  Offset  From  the  Plane  of  the  Resonators  by  a 

Distance  of  X/4. 


The  value  of  Gr  can  again  be  expressed  in  terms  of  ka  and  klno 
using  Equation  3-1,  where  Equation  3-3  is  used  to  express  klnm  in 
terms  of  n,  m,  and  klno  for  the  conical  configuration. 


klnm  =  y'r (klnof  -j|j"  2  sinj-J-  (n-mj)  where  n.m  =  1...N 


(3-3) 


The  gain  in  radiation  resistance  for  the  conical  configuration  for 
the  cases  N=4,  8  and  15  are  plotted  in  the  form  of  surface  and 
contour  plots,  as  shown  in  Figures  3-6a  through  6f.  Inspection  of  these 
plots  reveals  the  existence  of  the  same  local  maximum  as  the  planar 
configuration  of  magnitude  about  2  at  klno  =  n;,  and  ka  =  .01379.  Note 
that,  unlike  the  planar  case,  the  kl  axis  of  Figures  3-6a  through  f 
begins  at  kln0=1.75  vice  0.5.  This  is  due  to  the  displacement  of  the 
transducer  in  the  positive  z  direction.  This  displacement  requires  klno 
be  at  least  n/2. 
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The  local  maximum  values  for  Gr  at  klno  =  n,  ka  =  .01379  for 
various  values  of  N  ranging  from  2  to  25  and  their  associated 
parameters  are  tabulated  in  Table  3-2  for  the  conical  configuration.  In 
each  case  the  value  of  Gr  for  the  conical  configuration  is  slightly  less 
than  for  the  planar  configuration.  Also,  the  value  of  klmax  is  slightly 
greater  for  the  the  conical  configuration,  but  is  still  approximately 
equal  to  n.  (A  value  of  k  for  klno  corresponds  to  a  value  of  30  degrees 
for  the  angle  a  in  Figure  3-5.) 


N 

Gr 

k&max 

klnlax 

mag(UJf') 

arg(Uff') 

2 

1.27 

2.7 

.377 

80° 

4 

1.71 

1 

2.9 

.528 

78° 

6 

1.81 

3.1 

.421 

83° 

8 

1.82 

.01386 

3.1 

.321 

90° 

JO 

J  .82 

.01386 

3.1 

.255 

90° 

15 

1.81 

.01386 

3.r 

.168 

82° 

20 

1.82 

.01379 

3.1 

.128 

89° 

25 

1.82 

.01372 

3.1 

.103 

91° 

Table  3-2.  Values  of  Gr  for  Values  of  N  Ranging  from  N=2 
to  25  and  the  Parameters  Associated  with  the  Local  Maximum 
at  kl  s  k,  ka  =  klr  =  0.01379  for  the  Conical  Configuration. 
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d  ITS  3.3 

u‘  kl 

Figure  3-6c,d.  Three  Dimensional  and  Contour  Plots  of  Gr  as  a 
Function  of  ka  and  klno  for  N=  8,  Conical  Configuration,  1  atm. 
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Figure  3-6e,f.  Three  Dimensional  and  Contour  Plots  of  Gr  as  a 
Function  of  ka  and  klno  for  N=15,  Conical  Configuration,  1  atm. 


a  CONICAL  CONFIGURATION  BEAM  PATTERN 

To  demonstrate  the  improved  directivity  of  the  conical 
configuration  over  that  of  the  planar  system,  polar  plots  can  be 
de\cloped  using  the  axis  system  as  defined  in  Figure  3-7,  where  th~ 
transducer  is  displaced  is  in  the  positive  z  direction  from  the  plane  of 
the  resonators. 


Figure  3-7.  Axis  System  for  Beam  Pattern  Plots,  Conical 
Configuration. 

Figure  3-8a  shows  the  beam  pattern  of  the  conical  configuration  in 
both  the  y.z  and  the  z,0=45°  plane.  As  expected,  the  symmetry  in  the 
y,z  plane  about  the  z  axis  is  lost,  with  the  sound  pressure  level  (SPL) 
in  the  positive  z  direction  being  8  dB  greater  than  in  the  negative  z 
direction.  If  one  compares  the  beam  pattern  of  the  conical 
configuration  to  that  of  the  planar  configuration,  shown  in  Figure  3-8b, 
one  observes  that  the  SPL  in  the  positive  z  direction  for  the  conical 
configuration  is  1  dB  greater  than  for  the  planar  configuration  and 
lOdB  greater  than  for  one  transducer  alone.  Figure  3 -8b  clearly  shows 
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that  significant  directivity  can  be  achieved  with  only  a  minor  decrease 
in  total  power  gain,  Gr. 


z  axis 


Figure  3-8a.  Beam  Pattern  for  the  Conical  Configuration 
in  the  y.z  Plane  and  z,6=45°  Plane,  P0=  latm.  Scale  in 
Decibels  (dB). 
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Figure  3-8b.  Beam  Pattern  for  the  Planar  and 
Conical  Configuration  in  the  y,z  Plane,  P0=  latm. 
Scale  in  Decibels  (dB). 


3.  CONFIGURATION  DIMENSIONS 


The  results  of  all  calculations  performed  have  been  presented  in 
terms  of  the  dimensionless  radius  and  spacing,  ka  and  klno-  The  actual 
dimensions  for  the  radius  a  and  the  separation  distance  lmax  for  a  given 
frequency  are  obtained  using  Equation  3-4. 


a  =  ^L£L 

2n  f 


InO  — 


klnQC 

271  f 


(3-4) 


Taking  ka  =  kar  =  0.01379  and  klno  =  *  for  the  optimum  values  of  a 
typical  system  at  an  ambient  pressure  of  one  atmosphere  in  Equation 
3-4,  a  plot  of  the  resulting  dimensions  a  and  lno  as  a  function  of 
frequency  can  be  made,  as  shown  in  Figures  3-9a  and  3-9b. 
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Figure  3-9a.  Optimum  Resonator  Diameter  (a)  at 
Resonance  Versus  Frequency  for  a  System  at  P0=  lai 


Figure  3-9b.  Optimum  Resonator -Transducer  Separation 
Uno)  Versus  Frequency  for  Maximum  Power  Gain. 
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B.  ANALYSIS  WITH  CONSTANT  TRANSDUCER  ka,  AT  50 

ATMOSPHERES 

Two  major  assumptions  that  have  been  made  throughout  the 
analysis  up  this  point  need  to  now  be  modified.  The  first  is  that  the 
ambient  pressure  is  one  atmosphere  and  the  second  is  that  the  radius 
of  the  transducer  equals  that  of  the  resonators  (kao  =  kan).  In  the 
following  analysis,  an  ambient  pressure  of  50  atmospheres 
(corresponding  to  a  depth  of  500  meters)  will  be  taken  as  a  more 
realistic  value  for  an  operational  system.  Since  the  internal  acoustic 
impedance  of  a  bubble  is  a  function  of  depth,  the  value  for  kar  of  a 
single  bubble  at  resonance  will  shift  from  0.01379  (P0  =  latm)  to 
0.09754  (P0  =  50atm).  Also,  the  value  of  kao  will  no  longer  be 

considered  equal  to  that  of  the  resonators,  but  will  be  held  constant  at 
various  specific  values. 

In  the  previous  sections  it  was  shown  that  little  is  gained  when  the 
number  of  sympathetic  resonators  is  greater  than  about  four;  hence 
the  following  analysis  will  be  restricted  to  the  case  N=4.  Conducting 
similar  calculations  as  before,  contour  plots  and  surface  plots  have 
been  generated  for  the  cases  N  =  4,  kao  =  0.05,  0.10  and  0.20,  for 
both  the  conical  configuration  (Figures  3- 10a  through  3-1  Of)  and  the 
planar  configuration  (Figures  3-1  la  through  3-1  If).  As  in  the  previous 
analyses  where  the  ambient  pressure  was  taken  to  be  one  atmosphere, 
all  cases  display  a  local  maximum  of  GR  of  almost  2  at  around  ka  =  kar  = 
0.09754  and  kl  =  n.  These  maximum  values  are  tabulated  in  Table  3-3. 
Also  tabulated,  though  not  plotted,  are  the  cases  ka0  =  0.025  and 
0.150.  The  first  entry  in  table  3-3  for  either  configuration  was 


38 


calculated  holding  kao  equal  to  kan,  as  in  previous  calculations,  vice 
holding  ka0  constant.  A  review  of  Table  3-3  and  Figures  3-11  and  12 
demonstrates  that  the  value  of  kao  has  little  effect  on  the  achievable 
power  gain. 


configuration 

kao 

Gr 

ka^ax 

klmax 

BuEMaiial 

argCUlf1) 

planar 

kan 

BP 

.09998 

2.8 

0.657 

90° 

.025 

mm 

.09998 

2.8 

0.041 

90° 

.050 

1.88 

.09998 

2.8 

0.165 

90° 

.100 

1.89 

.09998 

2.8 

0.657 

o 

o 

.150 

1.89 

.09998 

2.8 

1.469 

90° 

.200 

1.89 

.09998 

2.8 

2.590 

90° 

conical 

ka„ 

1.78 

.09901 

3.1 

0.607 

95° 

.025 

1.78 

.09950 

3.1 

0.038 

95° 

.050 

1.78 

.09950 

3.1 

0.154 

95° 

.100 

1.78 

.09950 

3.1 

0.613 

95° 

.150 

1.79 

.09950 

3.1 

1.376 

95° 

.200 

1.79 

.09950 

3.1 

2.425 

95° 

Table  3-3.  Values  of  Gr  for  Various  Values  of  kao  Ranging  from 
.025  to  .2  and  the  Parameters  Associated  with  the  Local 
Maximum  at  kl  =  tc,  ka  =  kar  =  0  0.09754  for  both  the  Conical 
and  Planar  Configurations  at  P0=  50atm. 
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b.s 

Figure  3-10a,b.  Three  Dimensional  and  Contour  Plots  of  GRas  a 
Function  of  ka  and  klno  for  N=  4,  Planar  Configuration,  50  atm. 
kao  =  0.05. 
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Figure  3-10c,d.  Three  Dimensional  and  Co: 
Function  of  ka  and  klno  for  N=  4.  Planar  Co 
kao  =  0. 10. 
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Figure  3-10e,f.  Three  Dimensional  and  Contour  Plots  of  GRas  a 
Function  of  ka  and  kln0  for  N=  4,  Planar  Configuration,  50  atm. 
kao  =  0.20. 


Figure  3-lle,f.  Three  Dimensional  and  Contour  Plots  of  GRas  a 
Function  of  ka  and  klno  for  N=15,  Conical  Configuration,  50  atm. 


1.  BEAM  PATTERNS  AT  50  ATMOSPHERES 

Since  there  is  no  significant  variation  in  gain  with  the  value  of  kao, 
beam  patterns  were  only  developed  for  the  case  ka0=kar,  shown  in 
Figures  3-12  through  3-14.  These  plots  show  the  maximum  SPL  in  the 
positive  z  direction  for  the  conical  configuration  to  be  1.5  dB  greater 
than  for  the  planar  case  and  10.5  dB  greater  than  for  the  transducer  in 
absence  of  the  resonators. 

2.  CONFIGURATION  DIMENSIONS,  FIFTY  ATMOSPHERES 

The  actual  dimensions  of  the  resonator  radius  a  and  the 
separation  distances  lno,  for  a  system  operated  at  an  ambient  pressure 
of  50  atmospheres  were  calculated  using  Equation  3-4  Using  the 
optimum  value  ka  =  kar  =  0.09754  and  kln0  =  rc.  Since  the  values  of  ln0 
are  not  depth  dependent,  Figure  3-9b  still  applies.  Figure  3-15 
displays  the  relation  between  optimum  resonator  radius  and  frequency 
at  50  atmospheres. 
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Figure  3-12.  Beam  Pattern  for  Planar  and  Conical  Configuration 
in  the  y.z  Plane,  P0=  50atm.  Scale  in  Decibels  (dB). 

z  axis 


Figure  3-13.  Beam  Pattern  for  Planar  Configuration  in  the  y.z 
Plane  and  z,0=45°  Plane,  P0=  50atm.  Scale  in  Decibels  (dB) 
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Figure  3-13.  Beam  Pattern  for  Conical 
Configuration  in  the  y,z  Plane.and  z,0=45° 
Plane,  P0=  50atm.  Scale  in  Decibels  (dB). 
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Fig.  3-15.  Optimum  Resonator  Diameter  (a)  at  Resonance 
Verses  Frequency  for  a  System  at  P0=  50atm. 
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IV.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SUMMARY 

» 

The  problem  of  multiple  scattering  of  an  incident  plane  wave  from 
a  system  of  compact  resonant  scatterers  was  solved  using  both 
Tolstoy's  method  and  network  analysis.  It  was  shown  that  these  two 
methods  of  analysis  are  equivalent  and  that  both  predict  the  existence 
of  "quasiresonance".  Network  analysis  was  then  applied  to  a  circular 
array  of  scatterers.  where  the  incident  plane  wave  was  replaced  by 
spherical  waves  radiated  from  an  active  transducer,  located  on  the  axis 
center.  The  acoustic  advantage  of  this  system  was  quantified  by  the 
gain  in  the  radiation  resistance  seen  by  the  transducer  in  the  presence 

a 

of  these  scatterers  compared  to  a  lone  transducer.  Two  configurations 
*  were  analyzed,  (1)  the  transducer  in  the  plane  of  the  resonators  and 

(2)  the  transducer  displaced  by  one-quarter  wavelength  from  the 
plane.  Both  configurations  showed  a  gain  in  radiation  resistance  of 
almost  a  factor  of  two  and  improved  directivity  over  that  of  a  single 
transducer. 

B.  CONCLUSIONS 

The  numerical  analysis  of  a  transducer  surrounded  by  a  circular 
array  of  sympathetic  resonators  indicates  that  a  gain  in  radiation 
resistance  of  about  two  is  obtainable  compared  to  a  lone  transducer. 
The  physical  dimensions  of  such  a  system  operating  at  an  ambient 
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pressure  of  50  atmospheres  are  displayed  in  Table  4-1  for  several 
frequencies  of  interest. 


Frequency 

(Hz) 

Resonator  Radius  (a) 
(m) 

System  Radius  (R) 
(m) 

conical  planer 

10 

0.329 

64.9 

74.9 

100 

0.0329 

6.49 

7.49 

1000 

0.00329 

0.649 

0.749 

Table  4- 1 .  Dimensions  of  Planer  and  Conical  Configurations 
for  Various  Frequencies,  P0  =  50  atm.,  R  is  the  radius  of  the 
array. 

Although  the  calculations  presented  here  were  for  spherical 
bubbles,  the  results  also  apply  to  more  realistic  devices  (e.g.thin 
shelled  cavities,  balloons,  etc)  as  long  as  they  are  compact  (ka«l). 
M.  Strasberg  [Ref.  7]  showed  that  the  shape  of  a  resonant  body  has 
little  effect  on  its  fundamental  resonance  frequency.  For  example,  he 
calculated  that  an  oblate  spheroid  with  a  ratio  of  lengths  of  major  to 
minor  axes  of  four,  the  fundamental  resonance  frequency  differs  by 
only  a  factor  of  1.08  from  that  of  a  sphere  with  the  same  volume. 
Therefore,  any  deformation  of  the  resonators  from  a  true  spherical 
shape,  for  example  due  to  forces  such  as  buoyancy,  should  have  little 
affect  on  system  performance. 


C.  RECOMMENDATIONS 

The  results  of  the  research  reported  in  this  thesis  are  very 
encouraging.  It  is  recommended  that  follow-on  work  include  the 


design,  construction,  and  testing  of  an  array  to  verify  the  theoretical 
predictions  presented. 

Theoretical  analysis  remains  to  be  conducted  in  at  least  four  areas. 
First,  in  all  of  the  calculations  performed,  the  sympathetic  resonators 
were  assumed  to  be  gas-filled  air  bubbles  executing  radial  oscillations 
only.  A  more  realistic  model  of  a  sympathetic  resonator  may  be 
required,  based  upon  the  properties  of  a  practical  device.  It  is 
recommended  this  be  investigated.  Second,  Tolstoy  has  shown  that 
the  scattering  of  an  incident  plane  wave  from  a  linear  array  of  resonant 
scatterers  exhibits  strong  directionality.  Perhaps  a  linear  array  of 
sympathetic  resonators  offers  better  performance  gain  than  the 
circular  arrays  considered  in  this  research.  It  is  recommended  that 
such  a  configuration  be  investigated.  Third,  it  was  shown  that  in  the 
optimum  conical  configuration  the  ring  of  sympathetic  resonators  lies 
on  a  cone  with  a  apex  angle  of  approximately  120  degrees, 
independent  of  frequency.  This  suggests  that  a  broad  band,  constant 
directivity  system  could  be  constructed  using  a  series  of  concentric 
rings  of  resonators,  each  tuned  to  a  slightly  different  frequency  and 
placed  in  its  appropriate  location  along  the  surface  of  such  a  cone,  as 
indicated  in  Figure  4-1.  It  is  recommended  that  the  properties  of 
such  a  system  be  investigated.  Fourth,  in  the  calculations  presented, 
the  sympathetic  resonators  were  all  assumed  to  be  identical,  and  their 
locations  were  assumed  to  be  precisely  known,  as  well  as  that  of  the 
transducer.  It  is  recommended  that  the  effects  of  resonator 
nonuniformity  and  imprecise  location  be  investigated. 
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Transducer 


Resonator  r=a 


Figure  4-1.  Conical  Geometric  Configuration  Allowing  an  Expanded 
Frequency  Response  Using  Two  Concentric  Rings  of  Resonators  of 
Radius  a  and  a’,  at  a  Separation  Distance  from  the  Transducer  of  1  and 


APPENDIX 


PROGRAMING  INFORMATION 

Inquires  into  the  computer  programs  used  to  develop  the  system 
parameters  and  associated  graphics  should  be  addressed  to: 


Professor  S.R.  Baker 
Physics  Department,  Code  BA 
Naval  Postgraduate  School 
Monterey,  California  93943 
Tel.  (408)  646-2732/2729 
Autovon  878-2732/2729 
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